Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  449 1665 (N25E58 L=148)   418     57        0 ( 430/1639- 468/1691)  90.0 E59-E55 N23-N27
      2 1551 1605 (N42W29 L=235)   374     38        1 (1529/1581-1574/1626) -18.4 W26-W31 N39-N43
      3  337 1557 (N17E62 L=144)   435     60        1 ( 317/1531- 357/1582)  90.0 E64-E58 N14-N18
      4 1428 1425 (N25W19 L=224)   260     75        0 (1393/1400-1464/1446) -11.3 W16-W21 N24-N26
      5 1633 1332 (N24W35 L=241)   440    142        0 (1582/1273-1677/1379) -32.0 W30-W38 N20-N25
      6  139 1108 (S12E69 L=136)  3188    186        9 (  87/1045- 194/1192)  31.6 E80-E62 S16-S07
      7  323  806 (S30E48 L=158)  1266    195        3 ( 231/ 773- 421/ 835)   8.5 E58-E38 S32-S27
      8 1834  671 (S09W68 L=274)   626     61        0 (1800/ 643-1866/ 692)  29.7 W65-W72 S11-S07
      9 1470  507 (S28W43 L=248)  2666    509        4 (1269/ 448-1718/ 579)  31.0 W28-W60 S35-S17
  

    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