Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1524 1732 (N35W60 L=166)   826    157        1 (1462/1711-1573/1744) -32.0 W55-W63 N32-N37
      2 1034 1570 (N40W16 L=121)  1014    252        5 ( 980/1492-1095/1621) -84.8 W12-W18 N33-N44
      3 1468 1300 (N12W34 L=139)  2826    677        5 (1160/1081-1693/1462) -54.0 W19-W45 S05-N29
      4  454 1283 (N34E34 L= 71)   340     99        1 ( 422/1253- 491/1317) -38.7 E36-E31 N31-N36
      5  906 1123 (N15E05 L=101)   217     89        0 ( 876/1089- 940/1174)  33.7 E07-E01 N13-N17
      6 1646  863 (S17W36 L=141)   752    150        3 (1587/ 817-1684/ 902) -45.0 W31-W38 S20-S14
      7  530  512 (S14E44 L= 61)  1362    307        1 ( 504/ 396- 566/ 615)  47.3 E51-E38 S22-S08
      8 1357  449 (S37W09 L=114)  1248    235        5 (1309/ 389-1407/ 521)  74.7 W05-W12 S42-S31
      9  663  388 (S25E40 L= 65)   319     68        0 ( 628/ 371- 694/ 404)   0.0 E43-E37 S25-S23
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y